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Abstract 

We investigate the distribution of roots of polynomials of high degree with random coefficients 
which, among others, appear naturally in the context of "quantum chaotic dynamics". It is shown 
that under quite general conditions their roots tend to concentrate near the unit circle in the 
complex plane. In order to further increase this tendency, we study in detail the particular case of 
self-inversive random polynomials and show that for them a finite portion of all roots lies exactly on 
the unit circle. Correlation functions of these roots are also computed analytically, and compared 
to the correlations of eigenvalues of random matrices. The problem of ergodicity of chaotic wave- 
functions is also considered. For that purpose we introduce a family of random polynomials whose 
roots spread uniformly over phase space. While these results are consistent with random matrix 
theory predictions, they provide a new and different insight into the problem of quantum ergodicity. 
Special attention is devoted all over the paper to the role of symmetries in the distribution of roots 
of random polynomials. 
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1 Introduction 



Semiclassical approximations for multidimensional quantum systems and the manifestations of chaotic 
behaviour in quantum mechanics have attracted wide attention during the last years (see e.g. [|^, |||] 
and references therein). In these approximations one quite often needs to locate the roots of poly- 
nomials of high degree whose coefficients are rapidly-varying erratic functions of the energy. As a 
consequence, these coefficients may be considered as random variables, even in a small energy inter- 
val. Though the distribution of roots of polynomials with random coefficients have been studied in 
the past (see e.g. ||3|-|5|), its relevance with respect to the specific problems which naturally arise 
in the context of quantum chaotic dynamics (and in other domains of physics as well) has been 
underestimated, and a number of elementary and basic questions have not yet been solved. 

The main purpose of this paper is to study properties of random polynomials with emphasis on 
the above mentioned connection. A short version containing some of our results was already published 
in 1^]. Except for some general symmetry considerations, the coefficients of the polynomials will be 
considered independent random variables. This assumption is in some cases justified with physical 
arguments, and in others just made because of mathematical simplicity. Our investigations could 
have direct application in other fields. In fact, zeros in the complex plane of polynomials with random 
coefficients occur in a variety of problems of science and engineering (zeros of the partition function 
in statistical mechanics, theory of noise, etc). The present investigation may also be of interest if one 
views zeros of polynomials as interacting particles in two dimensions, as, for instance, eigenvalues of 
random asymmetric matrices can be physically interpreted as a two-dimensional electron gas confined 
in a disk @,§. 

The use of a statistical approach in the description of complex systems is an old idea. In particular 
the random matrix theory (RMT), originally formulated in the context of nuclear physics, has had a 
great success and impact in the study of quantum chaotic dynamic^l^ and disordered systems ||l|, ^J. 
As it will become clearer in the following sections, the statistical analysis of these systems by random 
polynomials is, in some sense, complementary to the RMT. 

In Section 2 we consider general random polynomials of degree N whose coefficients are independent 
random variables having zero mean. We show that under quite general conditions their roots tend to 
concentrate in an annulus near the unit circle of the complex plane, and that the width of this annulus 
goes to zero as ^ cxd. In |3| this result was proved by a different method for the particular case 
when all second moments of the coefficients are equal. Our method, based on the existence of a saddle 
point configuration, seems to be more general and physically transparent. 

^we will use this expression to avoid the more proper but lengthy term 'quantum mechanics of a classically chaotic 
system'. 
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Aside from section 2, which constitutes a general introduction, the paper can be divided into two 
main parts (sections 3 and 4) which are to a large extent independent. In Section 3 we investigate the 
important case of self-inversive polynomials, namely polynomials whose coefficients satisfy Eq. (|3.6| ) 
below and whose roots are hence distributed symmetrically with respect to the unit circle. This type 
of polynomials appear when considering either the semiclassical quantization via the transfer-operator 



method or quantum maps [11, We prove that for them a finite fraction of the roots lies 

exactly on the unit circle. Moreover, we compute analytically the two-point correlation function of 
these roots. We find that there is linear repulsion between them at short distances, and compare to 



the Gaussian Orthogonal Ensemble of RMT |13| , 14 1. The possibility of locating in a statistical sense 
all the roots on the unit circle is also discussed. 

In Section 4 we consider the properties of eigenfunctions of chaotic systems and their relation with 
the classical concept of ergodicity, in the spirit of Refs |l^, [l^. Using a phase-space representation, 
we write down the eigenfunctions of spin systems in a polynomial form and consider the coefficients to 
be random. We then analyse how the roots of these polynomials are distributed in phase space. The 
second moments of the coefficients now grow with N so fast that the resulting distribution of roots 
turns out to be uniform over the phase space, a sort of quantum ergodicity. We moreover compare 
this result with the predictions concerning the coefficients of the eigenfunctions of a random matrix 
ensemble. 

It turns out that the existence of a symmetry of the roots with respect to a line increases consid- 
erably the probability of finding roots exactly on that line. In subsection 4.2 we consider the influence 
of symmetries on the distribution of the roots of chaotic eigenstates. For that purpose we analyze 
the eigenfunctions of a certain quantum map having two antiunitary symmetries and show numeri- 
cally that their roots tend to concentrate over the associated phase-space symmetry lines. We also 
investigate how this phenomenon disappears as the symmetries are broken. 

In the appendices we explain in detail our computations. 

2 Some General Properties of Random Polynomials 

Given a distribution function 

P(ao, ai, . . . , aN)<f'aQ(f'ai . . . (fa^ (2-1) 



for some complex coefficients {ofc}, we are interested in the distribution in the complex plane of the 
roots {zfc}. A; = 1, . . . , A'^ of the polynomial 

P{z) = ao + aiz + a2Z^ + . . . + aNZ^ . (2.2) 
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For clarity, we consider the particular case where the real and imaginary part of the coefficients 
{cfc} are independent normally distributed real- valued random variables with zero mean and standard 
deviation ak (henceforth denoted as a GRI distribution) 

P(ao,...,a^) = n^exp(^-|^j . (2.3) 

The joint probability density for the zeros is obtained by changing from the variables {ao, oi, 03, . . . , cat} 
to the variables {ajv, Zi, Z2, ■ ■ ■ , zn} using the standard formulae 

QN-k = i-'^)'' aNUk{z), k = l,---,N 

^fcC-^) = Zi-^Zi2 • • • Zii^i 

l=il<i2<---<ik=N 

and integrating over ajv, with the result 

Dn{z ) 
[Gn{z)] 

where 



V{z^,...,zn) = Cn ,^ (2.4) 



GN{z) = l+Y.r]k\uk{z)\\ (2.5) 
k=l 

Dn{z) = n \zj - -ZfeP- 

j<k 

The factor D]\[{z) comes from the Jacobian of the transformation. 
Some properties of this distribution are: 

(i) if for all k ak = a , then 'D{zi, . . . , zn) is independent of a. This is a consequence of the fact 
that the roots of P{z) are unchanged by multiplying P{z) by a constant. 

(ii) T>{zi, . . . , zn) is invariant under a rotation of the coordinates in the complex plane, i.e. 

V{zi,...,ZN)=V{zie''^,...,ZNe''^) . 

(iii) if aN-k = cr^, then P(zi, . . . ,zn) is invariant under inversion with respect to the unit circle C 
in the complex plane, 

V{zi, ...,zn)= T^i^/zi, 1/zn), 
where the bar means complex conjugate. 
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In order to find the most probable distribution of roots, we have to locate the maximum of P. The 
equation for the extrema dlogD/dzp = can be transformed into the following form: 

N 

I mil \ 

QpUk = gp, p = l,...,N (2.6) 



k=l 



dz„ 



where 



dp 



One can easily check that the following configuration is a solution of these equations: 

2tt. 



Zk = rN exp 



1, 



(2.7) 



where 



rN 



N 



l/N 



JV^OD 1 / fJo 

- 1 + T7 In 

N \aN 



.N + 3, 

We will refer to this configuration as to the crystal solution. To prove that it obeys Eq.( |2.6| ) we remark 
that the given by Eq.(2.7) are solutions of the equation z^ = const, and therefore for this solution 



all Uk with A; = 1,2,...,A^ — 1 vanish and Eq.( |2.6[) reduces to 



VNUN 



dUN 

dz„ 



QpUN 



9p 



p=l, 



(21 



But duN / dzp = un/zp and gp = C/zp where 

N-l 



c 



1 



1 _ 1 (iV - 1) 

exp(i27rj7iV) ~ 2 (iV + 1)' 



Then Eq.(^) becomes 



rjN\uN\\l - C) = C ■ 



using the definitions of r/jv, un and C, the reader can verify that this equation is satisfied. Therefore 



all Eqs.(^^) for p = 1, . . . , N will be fulfilled by the crystal solution (2.7). 

It follows that (except in the case of an exponential (or faster) dependence of cjq/cjn on N) in 
the limit — > oo the radius r^r tends always to one. Moreover, the phase <j) in Eq.( |2.7| ) is arbitrary 
because of the property (ii). In the particular case ctq = ctn, tn is equal to one, as it must be according 
to property (iii). 

The existence of the crystal solution implies that the distribution of roots ^^{zi, . . . , zn) will always 
have an extremum on the line \z\ = r^r. But the actual distribution of roots in the complex plane 
depends, however, on the second derivatives of log P. There is a direct connection between their 
magnitude and the dominance of the crystal solution. For example, putting = gn and choosing 
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all the other cj's such that Nr]k ^0, [k = 1, . . . , — 1) as ^ oo, one concludes that the crystal 



solution will not dominate the distribution since the intermediate terms in the denominator of (2.4) 
automatically cancel, independently of the crystal solution. An extreme case of this behaviour is what 
we call, for reasons which will become clear later, the SU{2) polynomials, for which crk/(^N = \J C% 
where are Newton's binomial coefficients. As we shall see in section 4, for these polynomials the 
roots spread all over the complex plane. 

On the contrary, the case Nrj^ oo enforces the crystal solution since the factors in the 
denominator of (|2.4| ) are multiplied by an increasing function of A^. We shall consider in more detail 
several cases of this type in the next sections. 

These simple considerations show that under quite general conditions the roots of random polyno- 
mials of high degree for which the mean value of the coefficients is equal to zero will tend to concentrate 
near the circle \z\ = 1. In order to illustrate this with a numerical example, in Fig. la are plotted, 
in the complex plane, the roots of 200 trials of a polynomial of degree = 48 whose coefficients are 
GRI-distributed all having the same second moment. The observed distribution clearly satisfies the 
properties (ii) and (iii). 

On the other hand, it can be easily proved that if the mean value of the coefficients are non-zero, 
the roots tend to concentrate around the roots of the mean polynomial 



3 The Self-Inverse Symmetry 

As shown in the previous section and as illustrated in Fig. la, in the large- A^ limit and under certain 
circumstances the roots of random polynomials tend to concentrate around the unit circle but not, in 
general, on it. We now want to study some simple conditions to be imposed over the coefficients of 
a random polynomial in order to locate as much zeros as possible exactly on the unit circle. Among 
other areas of physics and mathematics, this problem is of interest in the context of semiclassical 
approximations in quantum mechanics because of the following reasons. 

Some recent methods to solve quantum problems incorporate the physical information into a 
transfer operator T which is an A^ x A^ unitary matrix |^] . The operator T has A^ complex eigenvalues 
{zk} , k = 1, . . . , N lying on the unit circle C of the complex plane, which are determined by the roots 
of the characteristic polynomial 

N 

P(z) = det(z - T) = ^ afc z'^ . (3.1) 

k=0 

The eigenvalues are functions of the energy of the system, = Zk{E), since T is a function of E. 
They move on C as -E is varied, and the quantization condition, which takes the form of a Fredholm 
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determinant, states that whenever one of the roots crosses the point z = 1 then the corresponding 
energy is an eigenvalue of the system 



det (1 - T{E)) = . (3.2) 
A similar equation but with no energy dependence is obtained when considering the quantization 



of maps, where T is now the one step unitary evolution operator of the map |11|. Since T is a 
finite matrix, it has N nontrivial invariants which can be chosen either as the coefficients {ofc} of the 
characteristic polynomial, or as the traces of the first powers of T 

N 

Ml = Tr T'^iE) = J2z^ , L = 0,...,N. (3.3) 

k=l 

One can express one set of invariants in terms of the other set by means of the following recurrent 
formulas 

1 

QN-k = -^"^aN-k+j Mj , k = l,...,N, (3.4) 

i=i 

with a^v = 1 by definition of P{z). These relations are of interest because in the semiclassical limit 
N ^ oo one can express the traces of the powers of T as an amplitude sum over the classical periodic 
orbits 

Here, ^{L) are all the periodic orbits corresponding to L iterations of the initial map at energy E 
(including repetitions), S.y{E) is the action of the periodic orbit, /i-y is the Maslov index and A.y{E) is a 
real function depending on the stability of the orbit. Inserting (|3.5| ) into (|3.4| ) we obtain a semiclassical 
approximation for each coefficient a^, now written in terms of sums over all the periodic orbits up 
to period k (the special combinations between them given by Eqs. (|3.4| ) are called pseudo-orbits). To 
compute all the coefficients, all the periodic orbits up to period are needed. 

Since the action S^{E) is a function of the energy, in the semiclassical limit the moments are rapidly 
varying functions of the energy and the coefficients ak are sums of products of these rapidly varying 
functions. It thus seems natural to adopt a statistical approach and to consider these coefficients 
as random variables. This way of proceeding differs from the usual statistical approach to complex 
systems in which, instead of the coefficients of its characteristic polynomial, the matrix elements of a 
relevant operator are assumed to be random. 

There is, however, an intrinsic difficulty when approximating the coefficients of the characteristic 
polynomial by Eqs.( |3.4|) and ( |3.5D , or more generally when considering them as independent variables. 
We are ignoring the correlations existing among them that guarantee the unitarity of T. (Without 
correlations and as was pointed out in section 2, the roots will lie close to C but not on it.) A 
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consequence of this is that the eigenvalues of T are no more located on C, and Eq.( |3.2D fails to 
determine the full spectrum (typically, some of the eigenvalues are missed). 

It may therefore be useful to find some simple conditions to be imposed on the coefficients {ak} 
in order to restore - at least partially - the unitarity of T. A necessary but not sufficient condition is 
the self-inversive (SI) property 

aN-k = exp (10) fflfc (3.6) 



which can easily be obtained from Eq.( |3.l[) factorizing the polynomial and making the substitution 

N 



Zk = exp (i^fc). is a real function of the energy, = vrA^ + J2 (^k{E)- In the semiclassical limit |10| 

k=l 



~ 7riV(S), 



where N{E) is the mean number of levels with energy less than E. Polynomials obeying Eq.(|3 
satisfy a functional equation 

P{l/z)=e-'^P{z)/z^ . (3.7) 
It follows from this equation that the symmetry (|3.6|) of the coefficients is reflected into a symmetry 



of the zeros: if z^ is a root, then 1/zk is also a root, i.e. the roots either lie on C or are symmetrically 
located under inversion with respect to it. 

From the semiclassical point of view, the advantages of imposing the self-inversive symmetry are 
twofold: firstly because it implements, in a simple manner, part of the unitarity of T; and secondly 
because it lowers the number of periodic orbits to be computed. In fact, we only need to know the 
periodic orbits up to the period N/2 instead of A^, since we only need to compute half of the coefficients 
(the others being determined by symmetry) (see e.g. |17, |l2|). 



3.1 The fraction of roots lying on C 

Being a necessary but not sufficient condition we don't know, however, how many zeros are located 
on the unit circle by the SI symmetry. In order to answer this question, we compute the fraction of 
zeros lying on C for SI polynomials of the form (remember that ajv = 1) 

N-l 

P{z) = 1+Y. ^kz'' + z^ , CN-k = ak , (3.8) 

k=l 

and consider the coefficients {a^} ,k = 1, . . . , {N — l)/2 as GRI distributed complex variables with 
arbitrary variances a^. For simplicity we consider the particular case = and consider N to be an 
odd integer, without loss of generality. By substitution z = exp(i^) in ( |3.8D , SI polynomials transform 
into real trigonometric polynomials: 

m = ^e-^^/^P (e^^) = cos (^e) + f: [ck cos [(f " A:) ^] + 4 sin [(f " A:) ^] } , (3-9) 
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where M = (N — l)/2, Cfc = 3f?e(ajt) and = Qm{ak)- The zeros of P{z) lying on C correspond now 
to the real zeros of the real function f{0). 

The average fraction of zeros lying on C is defined as 

1 /■^'^ 



< >-- 



N 



< pie) > de 



where < p{9) > is the average density of zeros on C 



p{9) = J26{e-0k) = 6[f{e)]\n0)\ 

k 



(3.10) 



(3.11) 



(Primes indicate derivative with respect to 9.) To compute < p{9) >, we use the method of Kac H] 
which exploits the following representations of 5[f] and \ f'\ 



(3.12) 



The advantage of such representations is that, when performing the ensemble average over the coeffi- 
cients 

< p{9) >= 'D{ao,ai, . . . ,aN) p{9) d^ao . . . d^ttN (3.13) 



the exponentiation of / and /' (who are linear function of the {o/c}) allows an easy computation of 
the integrals. 

The computation of < p(0) > and < > for arbitrary N and arbitrary second moments {(Jk} are 
straightforward but lengthy, and we include them in the appendix A. Here we give the result for the 
particular case of constant second moments = cr\/k and in the limit N ^ oo (see Eqs.( |3.20| ),( [3.28D 
below for the exact answer for arbitrary and arbitrary variances in the case of a SI polynomial of 
the form ( ^.19D ). We find that, to leading order in 1/A^, the average density of roots depends on a 
through the scaled parameter e = cr-v/iV, with the result 



N 

< P{^) >~ Y exp 



cos^ [^e] /e 



ire 



sm 



N 



+ 



vr 



^3 ^ 



dx exp 



-3sin^ ( 



Integration over 9 gives 



< lyie) >: 



1 



V2/e 



2ir J-V2/e 



dye 



TT 



dx / dip exp 



cos ip + 



3 sin^ ip 



x^ 



/{exf 
(3.14) 

. (3.15) 



To understand these results, consider first the limit e — > 0. In this case, the coefficients {a^} tend to 
have a very narrow distribution centered around zero, and the polynomial (|3.8| ) is well approximated 
by 

P{z) l + z^ ■ 
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the roots of this polynomial are given by Eg. ( p. 71) , with r^r = 1 and (p = tt/N. This is indeed the 
behaviour recovered from Eqs.(^]l4) and ( 3.15 ). In ( [3.15 ), when e — > the first term in the r.h.s. tends 
to one while the second tends to zero, implying < u >= 1. Moreover, in ( 3.1^ ) exp[— cos^(Y0)/e^] — > 
for any 9 except at 9 = ^ {k + ^) , k = 0, . . . , N — 1 where we get a delta peak, in agreement with the 
crystal distribution. 

On the other extreme, when e — > oo, we can neglect in ( ^.8]) the term 1 + and 

N-l 

P{z) ^ ~^ ^akz^ , ttN-k = ak ■ (3.16) 

fe=i 

In this limit, the first term in the r.h.s. of Eq.( p.l5| ) tends to zero, while the double integral of the 
second term tends to vr, and therefore 

< u{e) >'~°° I/VS, (3.17) 



while from ( 3.14 ) we recover a uniform density 

Eq. ( |3.17| ) answers the question of how efficient is the Sl-symmetry to locate roots of a random polyno- 
mial on C. We observe that indeed it has a strong effect in the distribution of roots, since it manages 
to locate a fraction l/\/3 ~ 57% of the roots exactly on the unit circle. (For a different proof of this 
result see Ref. [0].) 

Fig. lb shows the superposition of the roots of 200 iterations of a = 48 self-inversive polynomial 
with all the second moments equal. The total number of zeros is the same as in Fig. la; the strong 
concentration of roots on C is stressed in that figure by the reduction of the black intensity outside C. 

Fig. 2 displays the fraction of points lying on C as a function of e, Eq.( 3.15| ). For small e we observe 
the existence of a plateau which can be interpreted in the following way. At e = 0, as we said before 
the roots coincide with the crystal lattice and < >= 1. When e increases, and since the zeros are 
analytic functions of that parameter, they cannot immediately move outside C because this would 
violate the self-inversive symmetry (zeros come by symmetric pairs with respect to C). The only way 
they can get out from C is to first move along C until two roots become degenerate, and then split 
in the radial direction, one zero moving in the positive radial direction, the other towards the origin. 
The size of the plateau can be estimated as the typical perturbation needed to produce a coalescence 
of two roots starting from the crystal solution. 

It is instructive to compare our results to an analogous result due to M. Kac |^]. He considers the 
case of polynomials with real coefficients {a^} having a GRI distribution with all the second moments 
equal. These polynomials satisfy the functional equation 
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(the roots lie either on the real axis or come by symmetric pairs under reflection with respect to it). 
He computes the fraction of real roots and finds a much smaller effect of the symmetry as compared 
to Eq. dsTzl ), since he shows that < >~ In N/N as N ^ oo (he also considers the distribution of 
the zeros on the real axis, see ||5|). A numerical simulation of the distribution of roots of random 
polynomials with real coefficients is included in Fig.lc. The weak concentration of roots on the real 
axis can be appreciated in the figure from the fact that the density of points surrounding C is essentially 
unchanged as compared to Fig. la. For completeness, we plot in Fig. Id the distribution of roots of 200 
trials of = 48 Sl-polynomials with real coefficients. See also subsection 4.2 for an analogous result 
concerning SU(2) polynomials. 

To conclude this subsection, let us point out that there is a simple way to put all the roots of a 
self-inversive polynomial on the unit circle C (in a statistical sense). In appendix A we prove that the 
general formula for the fraction of roots lying on C for SI polynomials of the form 

N 

■P(^) = ^'^kz'' , aN-k = ak (3.19) 

fc=0 

where the {ak} are complex GRI distributed variables is given by 

2/52" .X 

<.>=-y^ (3.20) 

where 

N-l JV-1 

2 ^ / A^ \ ^ 

g,= Y,al , 92 = • (3.21) 

k=l k=l ^ 



This is an exact formula valid for arbitrary A^ and arbitrary variances cr| of the coefficients. The 
particular case of equal variances is explicitly written in Eq.( |3.28| ) below. This result was also obtained 



recently in [19|, where the reader can also find a geometrical interpretation of it. Consider now the 
parametrization 

where s is an arbitrary real number. Using the asymptotic expansion 

T.k^^^ + OiL'^), a>-l, 



from Eqs.( 3^2l| ) and ( |3.20| ) it follows, to leading order in A^ 



1 1 . X 

< >~ =, s> -- . (3.22) 

V s + 1 2s + 3 ' 2 ^ ^ 
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For s = (all the d's equal), we recover the previous result < u >= l/\/3. For s —>■ —1/2, then 
< I' >— > 1. For s = — 1/2 we can estimate the rate of convergence towards <i/>=lasA^— >oo from 



^k'^ ^ InN + C + 0{1/N) 

k=l 



(where C is Euler's constant). Then, from ( 3.2C| ) and ( |3.2l| ) we get 

3 1 



< z/ >~ 1 



+ 0{1/N), 



4/2. 



(3.23) 



(3.24) 



4 ln(iV/2) 

The convergence is therefore quite slow (for example, to put 98% of the zeros on C we need ~ 
3.8 X W^^). 



3.2 Correlations between the roots 

Having determined that the SI symmetry locates, in the large- limit, l/\/3 of the zeros on C, the 
next relevant question is what are the correlations existing between those roots. In particular, we 
would like to know if they repel each other or not and how their correlations compare to the random 
matrix theory (RMT). We have therefore computed, using the same techniques as before, the average 
two-point correlation function 

R2{r) =<p{e)p{9 + T) > 

for the set of roots lying on C for polynomials of the form ( 3.1S| ). For this kind of polynomials the 
function R2{t) does not depend on 6 since the distribution of roots is invariant under rotations. 

The exact result valid for arbitrary N and arbitrary variances is presented in appendix B (cf 
Eq.( [B.15| )). In the particular case of equal variances and in the limit ^ oo, r — > 0, Nt constant 



the two-point correlation function normalized to the square of the mean density ( 3.18| ) takes the form 

24 



M^) = R2i6)/iN/2TTV3f -- 
where 5 = tN/2ii and 

- - cos(7r(5) / 

■(^5)2 _ 



B arcsin f — ) + \/ 
A 



52 



(3.25) 



A = 


1 


f (^<5) 


8{7r<5)2 ■ 


3 


B = 


1 


cos( 




C = 


-I sin2(7r5) 



(7r(5)2 



sin(7n5) 1 sin(7n5) 



cos(7r5) 



sin(7r<5) 



siv? {-kS) 



(3.26) 



Fig. 3 is a plot of the function R2{5). For short distances there is a repulsion between the zeros. 
More precisely, from ( p. 25 ) it follows that 

R2[.o) ~ 



10^/3 
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This behaviour is reminiscent of the hnear repulsion between eigenvalues obtained for the orthogonal 
ensemble (GOE) of the random matrix theory, but with a different slope. However our computations 
correspond to systems without time reversal invariance, since the coefficients are complex. We therefore 
could have expected a quadratic repulsion, like in the GUE case. This particular point deserves further 
investigation. The long-range behavior shows on the other hand pronounced oscillations. 

We have also computed numerically the nearest-neighbour spacing distribution for the zeros lying 
on C, shown in Fig. 4. This distribution was investigated in Ref. [^], where a recursive procedure for 
p{5) was developed. Our numerical results are in agreement with those obtained in that reference. 

From the general results of Appendix B, it is also possible to compute the two-point correlation 
function for the case ak = k~^/'^ , when statistically 100% of the zeros lie on C. We now find that the 
correlations between zeros tend to be much closer to a crystal-like solution (i.e., strong oscillations 
that survive for large values of 5) than in the case of constant second moments. Fig. 5 displays the 
function R2{5) computed analytically from Eq.( [B.15| ) for a self-inversive polynomial with = k"^/"^ 
and N = 3601; for such a value of the fraction of roots lying on C is, from Eq.( |3.24 ), < v >= 0.9. 
In the light of the results of section 2, the emergence of a crystal-like behaviour when Uk = is 
to be expected since in this case Nrjk ^ as — > cxd. 



3.3 The variance of the number of zeros lying on C 

In section 3.1 we have established that SI polynomials of the form ( 3.19| ) have asymptotically on the 
average a fraction < u >= l/\/3 of their roots lying on C if the complex coefficients {a^} are GRI- 
distributed with the same second moment. Our purpose now is to compute the asymptotic behaviour 
of the variance of that average fraction, defined by 

fj^ =< 1/2 > _ < >2 . (3.27) 

From Eqs.( |3.20| ) and (B.16) the exact form of the average number of roots lying on C for arbitrary A^ 
and equal variances is 

^~ ~ 1 



<u >- 



and therefore 



On the other hand, 



^l/3 + l/Ar + 2/(3Ar2) 1 A ^ ^ ) + oH/N^^ 



< u - H . 

3 3N 



(3.28) 
(3.29) 



_ 1 
— TP 



' I^^'lt <p{e)p{e')>d9de' 



j^- < p{B) >d9 + j^- R,{9, e')de dd' 

^ + 2 /o2'^(2vr - r) i?2(r)dr] + 0(1/Ar2) 



(3.30) 



12 



where we have replaced the average number of roots on C by its asymptotic value N/ \/3 and we have 
also exploited the fact that R2{0, 0') depends only on the difference r between 9 and 9' to express the 
double integral as a simple one. Because R2{t) is symmetric with respect to r = vr, Eq. (|3.30| ) can be 
rewritten as 

2 1 47r /■'^ ^ , , , 

<Zy^>~^- + — ^/ R2{T)dT 



V^N Jo 

or, again normalizing R2{t) to the square of the asymptotic mean density R2{t) = i?2('r)/(iV/27r\/3)^ 



< >~ ^i— + ^ I R2{t) dr. 



(3.31) 



V^N 2,71 Jo 

As shown in Fig. 3, in the large- limit R2{t) has some oscillations on a scale r ~ 0{1/N) and tends 
to one for larger values of r. In order to separate the contribution to the integral in Eq.( ^.31 ) from 
the oscillatory part of R2{t), we rewrite it in the form 



R2{r)dT='^ r R2{5,N)d5 + 



2iTq/N 



R2{T,N)dT , 



(3.32) 



where 5 = tN/2'k and g is a parameter which is large compare to one but much smaller than N . In 
Eq.( 3.32| ) we have explicitly indicated the A^-dependence of i?2- Asymptotically (cf appendix B) 

R2{5,N) ^ R2{6) + 0{l/N) 
R2{t,N) ^ l + 2/Ar + 0(l/Ar2) 

where R2{S) is given by Eq.( |3.25 ). Therefore, keeping terms up to order 1/A^ 



2tt 



N 



2nq 27r 



27r 



R2{r)dT^-l i?2(5)d5 + vr-— 1 + - = ^ + - 1 



A^ A^ 



A^ 



l-R2{5) 



dS 



(3.33) 



Because i?2('^) tends to one as (5 ^ oo, we are allowed to take the limit (7 — > oo in the integral. 
Collecting Eqs.(g), ( |33lD , (^^) and ( pTj ) we finally get 



(^/3 



-A 



1 

iV' 



(3.34) 



with A = 1 — R2{S) d5. We were not able to compute analytically this integral, and made 
instead a numerical calculation. We obtained A ~ 0.44733. 



4 SU{2) Polynomials 

4.1 A theorem concerning the ergodicity of wave- functions 

In the previous section we were mainly concerned with the problem of locating on C as many as 
possible of the roots of a random polynomial and studying their correlations. This problem concerns 
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in particular the spectral statistics of asymptotic approximations of chaotic systems. In this section we 
explore a somewhat different and in some sense opposite problem: polynomials whose roots distribute 
uniformly (at least as — > oo) on certain surfaces. Our motivation is related to the asymptotic — > 
structure of quantum eigenstates of classically chaotic systems. 

Consider a system characterized by its total angular momentum J = {Jx, Jy, Jz) whose modulus 
J is conserved by the dynamics. The motion of the arrow J in the three-dimensional space can be 
represented by a point moving on the surface of a two-dimensional sphere, a Riemann sphere denoted 
S, which is in fact the phase space of the system. 

The equations of motion for J are such that the representative point is assumed to move on 5 in a 
chaotic way. This is possible if the angular momentum J is coupled to some external time-dependent 
field, typically a magnetic field. The simplest case is a periodic time-dependence, which we henceforth 
assume. Integrating the equations of motion of J over one period of the field, the classical dynamics 
for the point moving on the surface of the sphere reduces to a discrete map M acting on S 

j{n+l) ^ ^ l^jln)^ _ (4^^) 

These equations determine the position of the arrow at time t = n + 1 knowing its position at time 
t = n. 

The quantization of such a map introduces a one-period evolution operator U, the analog of the 
classical map M 

|^(n+l) >^ jj^^in) > ^ (4^2) 

where lip^"^^ > defines the quantum state of the system at time t = n. Because the modulus of J is 
conserved, then [U, J^] = and the Hilbert space is finite and {2J + l)-dimensional. We can choose as 
a basis of that space the eigenstates of Jz, Jz|?tt. >= hm\m >, m = — J, — J + 1, . . . , J. In particular, 
the eigenstates of the unitary operator U 

U\i;c.>=e'^"\i^a>, a = l,...,2J + l (4.3) 

can be written as 

J 

\i^a >= J2 > ■ 

m=—J 

The classical limit of such models corresponds to = 2J — > oo. For convenience we normalize 



the radius of the sphere, given by hy/ J{J + 1), to one. 

In order to have a unified semiclassical framework for both classical and quantum mechanics, it 
is convenient to introduce a phase-space representation of the Hilbert space. For that purpose, we 



project the eigenstates \ipa > into SU{2) coherent-states \z > [21, 22], ipaiz) =< z\ijja > with the 
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result (we drop from now on the subscript a) 

N 

ij{z) = Y,^akz'', (4.4) 

fc=0 

where N = 2 J, are the binomial coefficients and where we have shifted to the new label k = m + J . 



The complex variable z labeling the coherent states and appearing in the polynomial (iA) is 
connected to the variables (0,0) spanning the Riemann sphere by a stereographic projection of the 
plane onto the sphere, z = cot (0/2) e"*^. The function ip{z) is therefore an analytic function defined 
on the two-dimensional sphere. Because it is a polynomial of degree N, it has zeros in that space 
which completely determine (up to a global normalizing factor) the quantum state. 

Our purpose here is to analyze the structure of the eigenstates of C/ in a regime where the classical 
dynamics is dominated by chaotic trajectories; this means that classically the iteration of a typical ini- 
tial point covers in time the entire two-dimensional sphere 5 in a more or less uniform way. According 
to the correspondence principle, in the semiclassical limit the quantum states of such a system must 



tend, at least in a weak sense, to the microcanonical density [23|. The simplest asymptotic realization 
would be a function ip{z) whose modulus is uniform over S. This, however, is not an allowed solution 
since ip{z) has to have zeros in S, and therefore cannot be uniform. Moreover, the number of zeros 
proliferates in the semiclassical limit. The closest approximation to a uniform density would then 
be a function il^iz) whose zeros spread all over S. This behaviour was in fact already observed in 



Ref.[15| for eigenstates of chaotic systems. In the following, we make a precise statement concerning 
the ergodicity of the distribution of zeros for such systems. Assuming that the coefficients {a^} in 
Eg. ([4. 4]) are GRI-distributed, we prove that the zeros of ipiz) are indeed spread all over S and that 
their distribution is moreover uniform. 

Our assumption concerning the coefficients {a^} is motivated by the random matrix theory. As is 
well known, the statistical properties of the spectrum of classically chaotic systems are well described, 
in the universal regime, by the results of the RMT ||2^, Much less explored are the eigenstates of 
such systems, i.e. the statistical properties of the coefficients {a^} for chaotic systems and how they 
compare to the RMT (some results concerning this problem can be found in [^5|). The invariance under 
unitary transformations of the GUE ensemble of random matrices implies that the joint distribution 
function for the amplitudes must be 

TV 



\S- 



2{Af+l)l 



k=0 



where \Sn\ = 27r"/^/r(n/2) is the surface of a (n — l)-dimensional sphere of unit radius. When 
computing average properties of the zeros of ip{z), this distribution is strictly equivalent to a GRI 
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distribution (see ref. ||5|, p. 6) 



(2vr) 



N 



exp 




(4.5) 



Thus, the amphtudes {ofc} in RMT turn out to be gaussian uncorrelated random variables. 

Assuming this distribution function, our purpose now is to compute the associated average density 
of zeros on S, < p{0,ip) >, and compare it to the "ergodic" distribution conjectured above. Expressed 
in terms of the complex variable z, the density of zeros can be written (see appendix C) 



dz 



(4.6) 



The next step in the computation would then be to exponentiate both terms in the r.h.s. of Eq.( 
and compute the ensemble average over the coefficients. We find, however, that it is not necessary 
to exponentiate the Jacobian \dip/dz\^ . From Eqs.( |4.4| ) and ( [4.6| ), using the exponential expression of 
the delta function and computing the ensemble average we find, for arbitrary N, the result 

N d^z N 



< p{z) > d z 



TT (1 + |2|2)2 



47r 



sin 9 d6d(j) 



(4.7) 



i.e., a uniform distribution of the zeros on the Riemann sphere. The proof of Eq.(4/7) is given in 
appendix C. 

The result (^]^) constitutes a precise statement concerning the ergodicity of the (zeros of ) eigen- 
states of chaotic systems, and makes a connection between the concept of ergodicity and the RMT. 
Since the zeros are genuine wavefunction parameters, their equidistribution constitutes a stronger 
statement than just the limiting ergodicity of the Husimi function (which is a smooth quantity and 
bilinear in ip). 

As already mentioned, Eq.( [4.7D holds for arbitrary A^. In particular in the extreme case A^ = 1 
(J = 1/2) it means that the single root zi = — ao/ai of the monomial ^l^{z) = qq + a\z is uniformly 
distributed on the Riemann sphere if a\ and oq are complex variables having a GRI distribution and 
the same second moment. This result can be directly checked from Eq. (|2.3| ), 



P(ao, ai)d^ao(i^ai 
T>{z\^ a\)S z\S a\ 



Integrating over a\ we find 




d'^aad'^ai, 



ai\ d zid ai . 



vr(l + |zi|2)2 
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in agreement with (4^). A polynomial like (4^) generalizes this result to arbitrary N. The validity 
of the equidistribution of roots for arbitrary is however related to the gaussian nature of the 
distribution of the coefficients. We have numerically tested other distributions. For example, we find 
that for a uniform distribution of the coefficients in the interval [-1,1], the distribution of roots tends 
to be uniform only for large values of N. 

To conclude this subsection, let us mention that correlations between roots of random SU(2) 



polynomials were recently studied in |26] and compare to results obtained from quantum chaotic 



systems. Moreover, the general A;-point correlation functions were also recently computed analytically 



by J. Hannay |27|| . 



4.2 A physical example 

In order to illustrate the theorem (^]^) we consider a kicked-spin model, classically defined by the 
Hamiltonian 

oo 
n=— oo 

where /i and p are constant parameters. Integrating over a period At = 1 the equations of motion 
take the form of a discrete map 

J>+^) = R^ifi) (p4'A J». (4.8) 



Ri{X) represents a rotation around the i-th axis by an angle A, and J = {Jx, Jy, Jz)- The quantum 
mechanical analog of Eq. (|4.8| ) is the one-step unitary operator 

= e-^'^'^^e-^i'^' (4.9) 



acting on a (2 J -|- l)-dimensional Hilbert space. The stationary equation ([4.3| ) determines the eigen- 
phases coa and eigenstates {tpa > of U . According to the result Eq.( ^.7D , for parameters (/i,p) for which 
the classical map ( [4. 8] ) is dominated by chaotic trajectories we expect the zeros of the polynomials 
( [4.4D associated with the eigenstates of U to be uniformly distributed over the Riemann sphere. Before 
checking this, let us briefly mentioned some properties of U. 

The operator U is not generic since it has two symmetries: it commutes with two antiunitary 



operators |28| 



Ti = e^'^'^^e'^'^^^A:, T2 = e"'''"^"e'''-^^ A, 

where A is the usual antiunitary complex conjugation operator. These operators satisfy = = 1 
and the time reversal property TiUTi = T2UT2 = U^^. These two symmetries are nongeneric in the 
sense that (a) they are not just the conjugation operator usually connected to time-reversal invariance 
and (b) they depend on the parameter controlling, together with p, the dynamics of the system. 
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Because = 1, the existence of an operator Tj commuting with U imphes that the eigenstates of 
U can always be chosen T-invariant, T\ipa >= IV'a >• I^i the coherent-state representation, this latter 
equation imposes a functional equation on each polynomial ipai^) 



Mriiz))=Mz) i = l,2; a = 1, . . . , 2J + 1 , (4.10) 

where Ti, i = 1,2 are the classical versions of the quantum operators Tj which associates to each 
point of the classical phase space z an image Ti{z). Each of these transformations has a symmetry 
line, defined by the set of points z satisfying Ti{z) = z. In terms of the canonical conjugate variables 
{4>, cos 9) spanning the sphere, the two symmetry lines are given by 

sin 6 

cose = ± , ^ , (4.11) 

y ^ l±cos fi 

where the upper and lower sign holds for the Ti and T2 symmetry, respectively. 

Because of ([4.10| ), if Zk is a root of 4^{z), then Ti{zk) will also be a root. According to the results of 
section 3, if the roots are symmetric with respect to a line and if the coefficients of the polynomial are 
random, we expect a concentration of roots over that line. In the present case, we have two symmetry 
lines. Figure 6a shows the superposition of the 60 roots of the 61 eigenstates obtained by a numerical 
diagonalization of ( [4. 91) for J = 30, /i = 1 and p = Air, which classically looks fully chaotic [^]. We 
observe the expected concentration of roots over the two symmetry lines, a free-of-roots region close to 
them, and a tendency to cover in a more or less uniform way the remaining phase space (see ref. [^] for 
more details). Although we haven't computed analytically, we suspect that the concentration of roots 
over those lines is not macroscopic (i.e., asymptotically tends to zero), like the number of real roots 
of random polynomials having real coefficients. In this context, let us mention that it has been shown 
recently [^] that the asymptotic fraction of real roots of an SU(2) polynomial having real coefficients 
is 1 /\/iV. This should be compare to the log N/N fraction valid for the original problem proposed by 
Kac. 

In order to break both antiunitary symmetries Ti and T2 we add an extra term in the propagator 

U = e'^-^h-y^f^-^^e'^P'^^ . (4.12) 

Figure 6b and 6c show the superposition of the roots of the 61 eigenstates for the same p and fi as in 
Fig. 6a for t = 1 and t = 6, respectively. After a transition regime where the symmetry lines are still 
observed (even though the symmetry has been broken), for t = 6 the roots spread in a more or less 



uniform way over the whole phase space, as predicted in Eq.(4.7). 

Notice that although in order to have better statistics we have superimposed in the figures the 
roots of all the eigenstates, the ergodic theorem ( |4.7D holds for individual eigenstates. However, a 
direct numerical test based on a single eigenstate would need much higher values of J. 
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For other systems having different phase spaces, it is natural to expect - guided by semiclassical 



intuition - a result equivalent to Eq.(4/7), at least in the semiclassical limit. And indeed a related 



result was recently found [29| for the usual Bargmann representation |^0|] of quantum mechanics. This 
is a representation of Hilbert space in terms of entire functions 



Hz) = E 



k=0 



(4.13) 



The associated phase space is the two-dimensional plane labelled by the canonical variables {q,p), and 
z = (g — ip)/\/2. It was shown that if the coefficients {ofc} in Eq.( [4.13|) are GRI-distributed with all the 
second moments equal, the zeros of ip{z) are uniformly distributed over the whole plane. Moreover, 



if the sum in Eq.( 4.13 ) is truncated at a finite value A^, then the density of roots is uniform inside 
a circle of radius ^/N, and tends to zero outside. Generalizations to other geometries as well as the 



autocorrelation function of %Ij{z) have also been considered |29]. 



5 Concluding remarks 

We have studied the distribution and correlation of roots of random polynomials under several condi- 
tions. In section 3, we have explored, motivated by the problem of the semiclassical spectral properties 
of chaotic systems, different ways to increase the number of roots of a random polynomial located on 
C, and studied correlations between the roots. If all the coefficients of the characteristic polynomial 
have the same standard deviation, then asymptotically the self-inversive symmetry locates a fraction 
l/\/3 of the roots on C with a variance inversely proportional to the degree of the polynomial (section 
3.3). The two-point correlation function of these zeros behaves linearly for short distances (Fig. 3). 
However, if the standard deviations are not equal for all the coefficients but instead are given by 
Cfe = A;"-^/^, then the fraction of roots lying on C tends to one as ^ oo. But the convergence is slow 
(logarithmic), and the two-point correlation function much more crystalline-like (Fig. 5). Surprisingly, 
and unlike the case of standard polynomials with real coefficients (the problem of Kac) , we were able 
to compute explicitly the exact fraction of roots lying on C and their two-point correlation function 
for arbitrary A^ and arbitrary variances. 

To improve by less artificial means the number of roots lying on C we need to incorporate additional 
correlations between the coefficients. However no simple procedure exists, and the exact conditions 
for all the roots of the characteristic polynomial to lie on C are certain complicated determinantal 
inequalities for the coefficients (see e.g. Q]). 

The self-inversive property arises naturally in other problems of physics, as in the case of certain 
Ising models in statistical mechanics where the partition function takes a polynomial form when 
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written in terms of the fugacity |31]. There, the Sl-symmetry is connected to a spin- up spin-down (or 



particle- hole) symmetry of the system. As shown in |31], under certain assumptions all the necessary 
correlations implying unitarity are present and as a consequence all the roots of the partition function 
lie on C. In Ref. [^] the reader can find additional examples and references concerning the distribution 
of roots of partition functions in connection with the theory of phase transitions. 

Another famous example of concentration of the roots of a certain function on some simple curve 
in the complex plane is provided by the Riemann zeta- function C(^)- This function, which has no 
explicit random parameters entering its definition, satisfies the functional equation ,^(1 — z) = (,{z), 
where (,{z) = 7r~ ^ ^ '^T{z/ 2) (^{z). Accordingly, its roots are symmetric with respect to the critical line 
^e{z) = 1/2. The Riemann hypothesis asserts that all the nontrivial roots of (^{z) lie on that line. 

In section 4 we have shown that if we assume for the coefficients {ak} a GRI distribution with 
second moments equal to \Jc^ then the roots of the eigenstates of a classically chaotic spin system, 
Eq.( [4.4D , are uniformly distributed over the phase space (in this case, the two-dimensional sphere). 



20 



Appendix A 



In this appendix we compute the average density and average fraction of roots lying on C for a self- 
inversive polynomial of the form 

N-l 



P{z) = l+Yl ^kz"" + 



a-N-k = a-k 



(A.l) 



k=l 



for arbitrary N . The coefficients at are assumed to be complex independent variables having a 
Gaussian distribution 



P(ai, . . . jCm) 



1 / 1 \ 



(A.2) 



where M = {N — l)/2 [N is an odd integer). By the substitution z = exp(i^), the problem reduces to 
the computation of the average density and average number of roots of a real function in the interval 
< ^ < 27r (cf. Eq.Q) 

M 



f{9) = ie-'^^/2p(exp(i^)) = cos (^o) + 



k=l 



Cfc cos 



N 



k]e 



+ dk sin 



N 



k]e 



(A.3) 



where Ck = 5Re(afc) dk = 9m(afc). The derivative of this function with respect to 9 is 



no) 



N ,N 
— sin 
2 



(f.-|{(f-. 



Ck sm 



iV 



k 



N 



k I du cos 



N 



(A.4) 



The density of roots of f{9) is defined by 

p{9) = s[fm\f'm 

and using the Kac's representation ( p.l2 ) for 6[f] and |/'| we get 



Pie) 



27r2 



d^e 



dr] 2 

-oo V 



(A.5) 



In order to compute the average of p{9) over the ensemble ( |A.2| ), we proceed in the following way. 
We first replace f{9) and f'{9) in ( |A.5| ) by its definition Eqs.( [A.^ ) and ( |A.4 ). Then we average p{9) 
over the coefficients Ck and dk (it can be shown that the order of integration can be interchanged; see 
[kac]), 

/OO POO 
... d^ai . . . d^aMp{9)V{ai, ....au) ■ (A.6) 
-oo J — oo 

The final step is to integrate over rj and ^. 

1. We will need the average of e*^-^ e^"^^' . From (A.3) and (A.4) we have 



N 



N 



+ Vf' = Ccosi—9] - —r]smi—9] + {skCk + tfc4) , 



iV 



M 



(A.7) 



k=l 
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where 



Sk = ^ cos — kj ~ (t ~ ^)'? ^™ 
tk = Csin [(f - A;) + (f - A:)77cos 



(f - k) 
(f - 



From these equations and the fact that 

1 



2vra| 



(A. 



(A.9) 



we get from ([A.6| 



< e'^'^e'"'-^' >= exp <! i 



N 

icos ( —e 



N 



-r^sm 



— ( 
2 



1 f 1 



But si + tl = ^^ + {N/2 -kYri\ Then defining 

M 

91=11 (^l 



k=l 
M 

92= T, 

k=l 



N 



(A.IO) 



A; ai 



we can write 



< e^^^Q^'^f >= exp 



exp|i[cos(^^0) ^_^sin(^y0) r/]| . (A.ll) 



Note that the functions gi and 32 contain all the information concerning the variances of the random 
coefficients. For ?7 = 0, the result is 



< e^^-^ >= exp ( — (71^^ ) exp 



1 



N 

icos(-^)e 



Then from (|ATl| ), ( |AT^ ) and 



27r^ ' 



^fi_e-|W-if sin(fe)r, 
i-00 r?^ ^ 



(A.12) 



(A.13) 



2. The next step is to evaluate the integrals in ( [A.13| ). The integral over rj can be written 



m 



00 

-00 yf 



where /3 = (72/2. Since 



di_ 

d(3 



and, from ( |3.12| ) 



/(/3 = 0)= r^fl-e-'f-^T^)" 



vriV 



sm 



— ^ 
2 
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then 



m 



ttN 





r 








Jo 



Moreover, the integral over ^ in ( A.13| ) is 



/o ^/y 



2tt 



91 



(A.14) 



(A.15) 



Doing the change of variable x = \/y/P in the integral of Eq.( |A.14 ), the final result for the average 
density of zeros on C is 



< Pio) >-- 



-cos2(f e)/2gi f ^ 



sm 



N 



+ 2, 



3. The average fraction of zeros lying on C is defined as 

1 /■^'^ 



27r Jo 



^^g-iV2sin2(fe)/(8g2X^) 



1 r^^ 

<v>=-j^ < p{e) > de . 



(A.16) 



(A.17) 



The integral involving the first term between curly brackets in ( A.l(]| ) can be rewritten (by an obvious 
change of variables) 



2-K 



g -cos2(fe)/(25i) 



N 

sm I — i 



Ntt 



dip e 



■cos2(^)/{2gi)| 



sin^?! . 



The function to be integrated is periodic of period tt, and since sine/? is a positive function in that 
interval the previous equation takes the form 



2 / d(fe 2ffi '^"^ sincp = l^fgl 



dye 



-S/V2 



(A.18) 



The integral over the second term between curly brackets in ( |A.16 ) can be written, putting ip = NO/2 
and by the same argument as before 

„2,„ A72 „■ 2 



[\x /'"d0e-H'(f ^)/^^+^'^-'(f ^)/(4^2.2)]/2 = 2 /'dx rd^ exp 
Jo Jo Jo Jo 



COS If N sin tp 
91 4^2 a;2 



Putting together Eqs.( |A.16 )-( A.19), the fraction of zeros lying on C is given by 



< u > -- 



1 /-l/v^ 



dye y''/'^-\ — / dx I dipe 

vriV V 51 ^0 



2 fg^ 



2vr J-i/^ 



1 m 



2 I N sin ip 



(A.19) 
(A.20) 



Both results ( A.16| ) and ( A.20| ) are valid for arbitrary and arbitrary variances (t|. A; = 1, . . . , M. 
Moreover, both expressions are the sum of two terms: the first one in each expression is related to the 
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term 1 + z in ( |A.l| ), while the second comes from the random part of the polynomial. In the case 
where all the variances are equal, af, = cr^ Vfc, then the coefficients gi and 52 reduce to 



M 
k=l 

k=i ^ 



N 



(A.21) 



In the limit N ^ 00, gi (j'^N/2 = 52 — cr^A^^/24 = e^A^^/24, where we have introduced 



the scaled parameter e = ay N . Replacing these asymptotic expressions for gi and g2 in ( |A.16D and 



( A. 20 ), we find that < p{9) > and < z/ > depend only on the rescaled parameter e, and are given by 
Eqs.( 3.I4 ) and (|3.I5| ), respectively. 

On the other hand, when \ak\ » 1 VA;, then gi and g2 ^ 00 and we can ignore the first term in 



( A. 20 ). In this case (the pure random limit, which corresponds to neglect the term + 1 in P{z)), 
the fraction of roots lying on C, Eq.( A.20| ), reduces to Eqs.( 3.20 ) which in turn simplifies to Eq.( 3.17 ) 
if all the cr's are equal and if we keep only the leading term in Moreover, the average density of 

zeros tends to 

< p{9) >^ -.[^ . (A.22) 



Appendix B 



We compute here the average two-point correlation function 



i?2(r) =<p{e)p{e + T) > 



(B.l) 



for the roots lying on C of a self-inverse random polynomial of the form 



N 



P{z) = Y.akZ^ , 
fc=0 



where the coefficients , = 0, . . . , ^^^-2^ = M are complex independent variables having a Gaussian 
distribution ( |A.2| ) (we again assume for simplicity that N is odd; for even A^, we must take M = N/2). 



(B.2) 



(S.3a) 



Substituting z = exp{i6) in (|B.2| ) and ignoring prefactors we end up with the real function 

M 



/(^) = ^\ '^k COS 
A:=0 ^ 



A^ 



k]e 



+ dk sin 



N 



where = 5Re(afc), = Qm{ak), and whose derivative is 

M 



k=0 ^ ^ 



sm 



N 

— -k]dk cos 



2 



{B.3b) 
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The function f{d) has, as shown in appendix A, an average number A^/\/3 of zeros in the interval 
< ^ < 2tt. Moreover, since < p{6) > is in this case independent of 6 (cf. Eqs. (|A.22|) and (^Isl)), 
R2 depends on r but not on 9. Using the definition ( 3.11 ) of p{9) and Eqs. (|3.12| ), we can write the 
two-point correlation function as 



1 - e 



mf{e+r) 



(B.4) 



We proceed as in appendix A. We first replace Eqs.(B.3) in (B.4) and compute the average over 
the coefficients {ck,dk) according to the definition ( |A.6|) . Thenceforth we evaluate the integrals over 
and r]i. 

1. We will need < e'[«»^(^)+52^(^+^)+''i-^'(^)+''2^'(^+^)] >. From Eqs.(B.3) 

AI 

U{0) + 6/(0 + r)+ Vif'iO) + mf'iO + t) = J2 [(sfc + Uk)ck + (tk + Vk)dk] 

k=0 

where 

Sk = 6 cos [{N/2 -k)e]- {N/2 - k) r/i sin [(iV/2 - k) 9] 
tk = ii sin [{N/2 -k)9] + {N/2 - k)r]i cos [(iV/2 - k) 9] 
Uk = 6 cos [{N/2 -k){e + r)] - {N/2 - k) r/2 sin [{N/2 - k) {6 + r)] 
Vk = 6 sin [{N/2 -k){e + r)] + {N/2 - k) m cos [(A^/2 - k) {6 + r)] . 
Averaging over the Gaussian ensemble gives 



< e 



>= exp < --^[{sk + Ukf + {tk + Vkf]crl 
I k=0 



But 



{sk + Ukf + {tk + Vk? = + il + {N/2 - k?{r]f + r?i) 

+2 [6^2 + {N/2 - k fmV2] cos [{N/2 - k) r] 
+2(iV/2 - A:)(6r/i - iim) sm [(iV/2 - A:) r] , 
which is independent of 9, as it should be. Then defining 

M 

91= E (^l 

k=0 
M 

52 = E iN/2 - kfal 

k=0 
AI 

53 = E cos[(iV/2 - k)T]al 

k=0 

AI 

94 

k=0 
M 

32 „ /a_2 



d93/dT = E iN/2 - k) sin[(iV/2 - k)T]al 



55 = -d'gz/dr' = E {N/2 - kf cos[(iV/2 - k)T]al 



(B.5) 



(B.6) 



A;=0 
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we obtain the result 



< exp {i[U{0) + af{0 + t)+ VifiO) + mf'iO + r)]} >= 

exp {-[giiCl + Ci) + 52(r?f + r/i) + 2(5366 + g^mm) + 254 (6m - 6??2)]/2} 



(B.7) 



From this expression the average of the different terms appearing in (B.4) can be evaluated, with the 
result 



|l _ e-(92'??+2g46m)/2 _ g-(92»?i-2g46'?2)/2 _^ (^-[g2(vl+V^)+2g5ViV2+294((2Vi-^iV2)]/2'j ^ 



2. The next step is to compute the integrals in (B.8). The integrals over 6 and 6 are straightforward, 
since they involve exponentials of quadratic forms: 

1 



2^3,./.2_„2 



where we have introduced 



drji dr]2 
111 V2 



gigl 



2\ ' 



/3=<?5 



g3 gj 



Consider now the integral 

m = 

and write it in the form 



^ drji dri2 

00 f?i 



/3=0 



From Eq.( |A.15 ) evaluated at = we get the first term of the previous equation 

drix dri2 



m\ 



13=0 



00 Vi V2 ^ /V / 



Moreover 



dp 



dm dr]: 



'2 -^{vlW2)f,-l3viV2 



-e 2 > 



00 Vi V2 



implying dI/d(3\g^Q = by antisymmetry. Furthermore, 



27r 



x/a2 - /52 



and integrating twice 



( dx dy 
JO 



27r 



2-K 



(3 arcsin(/3/a) + \ja^ — ffl — a 



(B.9) 



(B.IO) 



(B.ll) 



(B.12) 



(B.13) 



(B.14) 
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We thus obtain, substituting in ( |B.12 ) the different terms 

a 



m = 



P arcsin 



Then the exact average two-point correlation function, vahd for arbitrary N and arbitrary variances, 
is 



R2(t) 



1 



/3 arcsin + \J — (3"^ 



(B.15) 



The sums 0.6D - needed in the computation of the coefficients a and (3 - can be evaluated explicitly 
if the variances are equal. We find 

gi=a\N + l)/2 
c,2 = {N^/2A + A^VS + iV/12) 
93 = cj^ sin [t{N + l)/2] / [2 sin(T/2)] 
54 = -T^) [n-os [r{N + l)/2] - 



(B.16) 



55 



( sin(TjV/2) cos(t^ 
4 \ sin(T/2) 



/2) 



Ar2 



1 



sm2(T/2) 



+ iVcos(riV/2) 



+ 



1 



sm^(T/2) 



In this latter case, taking the limit — > oo, r ^ 0, A^r — > constant, and normalizing R2 to the 
square of the asymptotic mean density (|3l8|) , Eqs.(pl5|), ( pl6| ) and (^l0|) can be rewritten in the 
form 



Appendix C 



We prove here the theorem Eq.( [4.7D , stating that the average density of roots of a polynomial of the 
form 



N 



(C.l) 



fc=0 



is uniform over the Riemann sphere for arbitrary N. In Eq.(C.l) the coefficients ak, k = , . . . , N are 
assumed to be complex independent variables having a gaussian distribution with the same standard 
deviation ak = <J \/k, and the are the binomial coefficients = N\/k\{N — k)\. 

We want to compute the density of zeros of a complex function in the complex plane. By definition 



p{z) = 5[^e ii^iz))] (5[9m (V'(^))] 



9SRe(V>) aSRe(V>) 
d^e{z) d'iTn{z) 



(C.2) 



By the Lemma 7.1, page 150 of Ref.[^ concerning the Jacobian of complex analytic functions, ( |C.2| ) 
can be rewritten 



p{z) = <5[5Re iipiz))] 5[Qm iij{z))] 



(hp 
dz 



(C.3) 
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We will use, for convenience, polar coordinates in the z-plane 



z = re 



lip 



and write = + idk- Then, from ( |C.1| ) 



N r I ^ (C.4) 

fc=o I 



Moreover 



dz 



N 



fc=o 



and then, in polar coordinates 

2 N 



dip 



dz 



kl=0 



(C.5) 



To compute the average over the Gaussian ensemble we will only exponentiate the delta functions in 
( p.3| ), but not the Jacobian. Using the expression ( |3.12| ) for the delta functions, and ( C.4D -( pl5| ) for 
/ and the Jacobian, respectively, the density (p.3|) can be expressed as 



N 



^ ' k=0 



N 



N 



+ E \/C^C^ki [ckCi + dkdi + iidkce - Ckdi)] rk+i-\Kk-i)p Uxp <; ^ (a„c„ + ) 

k=/=e=0 ^ J [n=0 J 

(C.6) 



where 



an = i^/C^r"'[cos{n^p)^l + sin(nv9)^2] 
= i^/C^r''[cos{nip)^2 - sin(n99)^i]. 
1. The average over the coefficients Ck and dk in (|C.6D involves expressions of the type 



(C.7) 



< cl^df exp < {anCn + I3ndn) \ > 
[n=0 J 



where the symbol < . > represents the average over the ensemble (2.3) taking all variances cr^ equal; the 
parameters ji in the latter expression can take the values 1 or 2. The computation is straightforward; 
for example, for ji = 2 and j2 = 

< 4 exp I E(«"'^« + /^"^«)f >= ^^(1 + «fc^^)exp I Y E("n + • 

[n=0 ) I n=0 J 



From ( |C.7D we can write 
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and hence 



N 



< 4 exp ^jy^^nCn + Pndn)j >= a\l + afa^) exp |- — (1 + r^iCf + C?) j . 
The other averages are computed analogously. The result of averaging (C.6) is 

9 r 9 1 { N r 

< p{r, ^)>= j^J JZo dCid^2 exp [-4(1 + r')^{e^ + Ci)] E C^^^' [2 - a'C^^r^Hei + Ci^ 



(C.8) 



2. The integrals involving the first term between curly brackets in ( C.8 ) give 



[2 - a'C'^r'Hei + Ci)] exp {-^^(1 + ifi + ei)/2} = ^ 



k ^2k 



(l+r2)^ 

Moreover, from Eqs. (|C.7| ) it follows that 

Using this result, we evaluate the integrals involving the second term between curly brackets 

4lT 



(l + r2)^ 
(C.9) 



Using Eqs.(|C]|) and ( |cl^ ), (Q can be expressed as 



cj4(l+r2)^ 



1 



7r(l +r2) 



2\N 



N 



k=0 



(1 +r 



2\N 



' N 



k u^2k~l 



(C.IO) 



(C.ll) 



\k=Q 



In order to evaluate the sums, we consider the identity 

N 

^N2k 



5:C^V- = (l + r 



2\N 



(C.12) 



fc=0 



By differentiating Eq.( p.l^ ) once and twice with respect to r"^ we get 



N 



^ C%k r^^-^ = r^^^^±^ =rNil + r^^^-^ 



k=0 



9(r2 



and 



fc=0 



(C.13) 



(C.14) 



r2A^(iV - 1)(1 + r2)^-2 + A^(l + r^) 



^2\N-l 
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respectively. Substitution of ( p.l3| ) and (p. 14 ) into ( |C.ll ) gives, finally 



< p{r, if) > rdrdif 



N 



N 



(C.15) 



vr (l+r2)2 ^ TT (l + |z|2)2 ■ 

A simpler form for the density of zeros is obtained projecting (by a stereographic projection from 
the north pole) the complex plane into the two-dimensional Riemann sphere (having unit radius), 
spanned by the spherical variables {9, (p) 



z = cot {e/2) e'^ . 



This transformation explicitly shows that, in fact, the density (C.15) is uniform on that surface 



N 

< p{9, (p) > dOdip = — sin 6 dOdip . 

4tt 



(C.16) 
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FIGURES 



FIG. 1: Distribution in the complex plane of the roots of random polynomials. In all the ((a) to (d)) 
cases we have superimposed the roots of 200 different trials of a polynomial of degree = 48 
whose coefficients obey a GRI distribution all having the same second moment a. (a) Standard 
random polynomial with all its coefficients complex and independent, (b) SI random polynomial 
with complex coefficients; only half of the coefficients are random, the other half being determined 
by complex conjugation, (c) standard random polynomial with real coefficients, (d) SI random 
polynomial with real coefficients. 

FIG. 2: The asymptotic average number of roots lying on the unit circle C for a SI random polynomial 
of the form (|3.8D as a function of the parameter e = a^/N . 



FIG. 3: The asymptotic two-point correlation function i?2 (Eg. ( p. 25] )) for the roots lying on C of SI 
polynomials with complex GRI-distributed coefficients. 

FIG. 4: Nearest-neighbour spacing distribution for the case of Fig. 3. 

FIG. 5: The two-point correlation function R2 for the roots lying on C of A'^ = 3601 SI polynomials 
with complex GRI-distributed coefficients and second moments = A;"^/^. 

FIG. 6: Phase-space distribution of the zeros of eigenstates of a chaotic system. In the three parts 
of the figure we have superimpose the 60 roots of the -|- 1 = 61 eigenstates obtained by 
numerical diagonalization of the kicked-top map ( 4.12| ) for ^ = 1 and p = Air. (a) t=0, where 



two antiunitary symmetries exist. We observe a concentration of roots over the two associated 



phase-space symmetry lines given by Eg. ( 4. 11 ). (b) t=l, both symmetries are now broken 



but there is still a concentration of roots on the symmetry lines, (c) t=6, all vestiges of the 
symmetries have disappeared, and we recover a uniform-like distribution, in agreement with the 
prediction (p^). 
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